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CMB Anisotropics from Outflows in Lyman Break Galaxies 



NO 
O 
O 
(N 

Q-i 

<D 
00 



> 
ON 
ON 

m 
On 
O 

NO 

o 

Oh 
i 

o 



13 



Daniel Babich 1 ' 2 * and Abraham Loeb 1 ! 

1 Astronomy Department, Harvard University, 60 Garden Street, Cambridge, MA 02138, USA 

2 Current address: California Institute of Technology, Theoretical Astrophysics, MC 130-33 Pasadena, CA 91125, USA 



5 February 2008 



ABSTRACT 

Thomson scattering of the Cosmic Microwave Background (CMB) on moving elec- 
trons in the outflows of Lyman Break Galaxies (LBGs) at redshifts 2-8 contributes to 
the small-scale CMB anisotropies. The net effect produced by each outflow depends 
on its level of deviation from spherical symmetry, caused either by an anisotropic en- 
ergy injection from the nuclear starburst or quasar activity, or by an inhomogeneous 
intergalactic environment. We find that for plausible outflow parameters consistent 
with spectroscopic observations of LBGs, the induced CMB anisotropies on arcminute 
scales reach up to ~ lfjK, comparable to the level produced during the epoch of 
reionization. 

Key words: cosmology - cosmic microwave background - galaxies :high-redshift 



1 INTRODUCTION 

Several experiments to observe the Cosmic Microwave 
Background (CMB) anisotropies on arc minute scales are 
currently, or wi l l soon be, un derway jKosowskvl l2003t 
iRuhl et all 12004 ILTet alJl2005l) . These experiments plan 
to measure the CMB power-spectrum for a spherical har- 
monic multipole index of 10 3 ^ I ^ 10 4 at several frequen- 
cies centred around the thermal Sunyaev-Zel'dovich null of 
217 GHz. Photon diffusion damps the CMB anisotropies on 
these small sc ales during cosmological recombination at red- 
shift z ~ 10 3 llSiUd ll968l. and so any observed signal must 
originate at much lower redshifts. Indeed, the above exper- 
iments plan to constrain the epoch of reionizat ion and the 
growth of structure in the low redshift universe Kahn et alJ 



flows in Lyman Break Galaxies (LBGs) to the small-scale 
CMB anisotropies. LBGs are believed to be the ancestors of 
present-day luminous elliptical galaxies. They are observed 
to produce gas outflows with velocities of several hundred 
km s _1 (see iGiavaliscol J2002) for a comprehensive review). 
In contrast with the traditional kinetic Sunyaev-Zel'dovich 
effect where the bulk velocity of the virialized gas is respon- 
sible for the induced CMB anisotropy, we focus here on the 
Doppler effect of the outflowing gas and ignore any bulk mo- 
tion of the LBG as a whole (which produced a smaller effect 
at the redshifts of interest). This bulk effect is included in 
standard calculations of the non- linear generalization of the 
Ostiker-Vishniac effect <Hull2000lh 

The contribution from a single LBG to the fractional 
temperature fluctuation of the CMB can be expressed as 



The primary physical mechanism which is responsible 
for the small scale CMB anisotropies is Thomson scatter- 
ing of CMB photons off moving electrons. Any peculiar 
velocity induces a Doppler anisotropy of the scattered ra- 
diation along the direction of motion. This accounts for 
the CMB anisotropies produced by the peculiar veloci- 
ties of clusters (the so-called kine tic Sunyaev Zel'dovich 
effect) JZel'dovich fc Sunvaevllliiol) . by peculiar velocities 
of linear overdensities in the intergala ctic medium (the so- 
called Ostriker - Vishniac effect) flOstriker fc Vishniadll98rl 
IVishniadlT987h and by the peculiar velocities of the fluctu- 
ations in the ionizatio n fraction during patchy reionization 
jGruzinov fc Hulll99Slh 

In this Letter we examine the contribution of out- 
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where <jt = 6.65 x 10 - 25 cm 2 is the Thomson cross section, 
n e is the electron number denisty, v is the electron pecu- 
liar velocity, c is the speed of light and h is the observer 's 
line-of-sight toward the LBG jZel'dovich fc Sunvaevlll98Clh 
The integration traces the photon's path through the LBG 
outflow. 

The radial extent of the outflow is found by solving 
the corresponding hydrodynamics equations. These coupled 
non-linear partial differential equatio ns can be reduced to 
a single ordinary differen tial equation jTegmark et alj[l993l : 
iFurlanetto fc Loebl2003r) under the assumption that the gas 
swept-up by the outgoing blast wave lies in a thin shell be- 
hind the propagating shock front (the so-called thin shell 
approximation). The validity of this approximation is illus- 
trated by the self-similar Sedov-Taylor-von Neumann solu- 
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tion for a point explosion in which 90% of the swept-up 
mass resides in a shell of thickness 10% of the outflow's ra- 
dius (see Ostriker & McKee 1988 and Ikeuchi et al. 1983 
for additional discussion on this approximation). The thin 
shell approximation allows us to treat the radiative transfer 
of CMB photons through the shock front in a plane parallel 
geometry. When the thickness of the shock front is small 
compared with the shock front's radius of curvature, the 
path length through the shock front can be expressed as 
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where 8R is the thickness of the shock front. Using conserva- 
tion of mass and the density compression ratio for a strong 
adiabatic shock, one gets 8R/R = (7 — l)/(37 + 3), where 
R is the r adius of the outflow and 7 is the adiabatic index 
of the gas JOstriker fc McKedll983) . 

Within the thin shell approximation the line-of-sight 
integration in Eq. Q is simplified to 
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[n e (ri)v(fi) - n e (r2)v(r 2 )], (3) 



where ri and r*2 are the location, in a coordinate system cen- 
tred on the LBG, on the shock front where the line-of-sight 
respectively enters and exits the shock front. The subscript 
i labels the contribution from a particular LBG. 

If the outflow is spherically symmetric then there is an 
exact cancellation of the anisotropies produced at the en- 
try and exit points of the line-of-sight through the LBG 
shock front 1 . However, observations of low-redshift star- 
burst galaxies, which serve as analogs of the higher red- 
shift LBGs, show evidenc e for highly non-spherical outflow 
geometries jMartinlll999l) . It is therefore reasonable to ex- 
pect that the early stages of LBG outflows, whether they 
are driven by starburst or quasar activity 2 , would produce 
CMB anisotropies. 

For simplicity, we will assume that the outflow is ax- 
isymmetric about some axis z (possibly the rotation axis of 
the galactic disc or possibly the jet axis of a central quasar) 
and expand both the outflow velocity and shock front den- 
sity in Legrendre polynomials, v(r) = f ~}2 e ve(r)Pi(z ■ r) 
and n e (r) = £\ ne(r)P((z ■ f), where we have assumed that 
the velocity field is radial (v = r). Even if the outflow began 
highly collimated, it will eventually isotropize as it propa- 
gates into the intergalactic medium (IGM; an analogous ten- 
dency exists in relativistic flows, see Ayal & Piran 2001). On 
longer timescales, IGM inhomogeneities, such as filaments 
and voids, will once again make the outflows non-spherical. 
Numerical simulations are needed in order to properly model 
these effects. Here we parameterize the level of asphericity 
in the outflow by a coefficient e (see Eq.Q below). 



1 Note that even if the outflow had a perfect spherical shape, 
the finite light crossing time through the outflow would pro- 
duce a net non-zero signal because the flow parameters are time- 
dependent. This effect would produce a signal a-^n e R{h.v) / c ~ 
<TTn e (dv/ dt)(R/ c) 2 that is extremely small and is ignored here. 

2 In this Letter we will only include the feedback driven by su- 
pernovae. Energy input from quasars will le ad to enhancements 
in the bubble size and outflow velocity [see iFurlanctto fc Loebl 
1200 ifl for a description of quasar outflows] . 



The outline of this Letter is as follows. In §2 we ana- 
lyze the CMB anisotropy induced by a single LBG with an 
arbitrarily non-spherical outflow. In §3 we calculate the re- 
sulting CMB power spectrum. In §4 we present numerical 
results and in §5 we summarize our conclusions. Through- 
out our discussion, we will a ssume the WMAP3 cosmological 
model JSpereel et aUl200db 



2 SINGLE LBG SIGNAL 

The signal along a given line-of-sight includes contributions 
from LBGs with different masses, outflow ages, orientations 
of the outflow symmetry axis and line-of-sight impact pa- 
rameters. For simplicity, we begin by analyzing the expecta- 
tion value for LBGs formed at a certain redshift and of cer- 
tain mass and age. Since the net effect from a single LBG can 
be either positive or negative we will find that the mean of 
the signal along a given direction is always zero, but a non- 
zero variance will be produced due to Poisson fluctuations 
in a manner equivalent to a random walk. The analysis will 
be done separately for the distinct cases where either one or 
two lines of sight intersect the same LBG outflow. 



2.1 One Sightline 

First we consider the case where one line-of-sight intersects 
a single LBG outflow. Averaging over symmetry axis orien- 
tation and impact parameter, the expectation value of the 
fractional temperature perturbation is 
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Here z is the symmetry axis of the outflow and b is the 
impact parameter at which the line-of-sight enters the LBG 
(sec Fig.0for definitions of the variables we use in describing 
the LBG outflow). We assume uniform probability distribu- 
tions for the symmetry axis orientation, P[z) = 1/47T, and 
for the impact parameter, P(b) = 1/irR 2 . 

Now we change the integration variable from the impact 
parameter b to the location of entrance point of the line-of- 
sight into the LBG outflow r. Thus, the expectation value 
defined in Eq. (3) can be written as 

A T J J h-K 4vr tv \ T J. 

where the integration over the location of the impact pa- 
rameter is restricted to a single hemisphere. Performing this 
integration we find that the signal vanishes on average, as 
expected from the fact that the net signal for a given LBG 
is just as likely to be negative as it is to be positive. 

There will be a non-zero contribution from a given LBG 
since the induced anisotropies will be different at the entry 
and exit points of the line-of-sight through the shock front. 
Because the LBG symmetry axes is randomly oriented with 
respect to the direction of the observer, the signal vanishes 
once the average over this direction is done. This implies 
that the mean signal and therefore the one-point function 
vanishes. We are ultimately interested in the two-point cor- 
relation function and the related power spe ctrum. There are 
two c ontributions to the power spectrum ijCoorav fc Shetbl 
2002). The first is a clustering (two halo) term, originating 
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Figure 1. Schematic illustration of the geometry for a single LBG 
outflow. The black circle represents the shock front, the orange 
arrow labelled n is the line-of-sight from the observer, the red 
arrows labelled ri , f*2 are the entry and exit points of the line-of- 
sight on the outflow, the blue arrow labelled is is the symmetry 
axis of the outflow and the green arrow labelled b is the impact 
parameter of the line-of-sight. 



from the correlated perturbations in the cold dark matter 
density produced during inflation. The second is a Poisson 
(one halo) term, originating from Poisson fluctuations in the 
number density of halos. The one-point function vanishes be- 
cause net temperature anisotropy produced by a given LBG 
is uncorrelated with the signal from other LBGs along the 
line-of-sight. The clustering term only implies that the num- 
ber density of LBGs nearby another LBG is greater than av- 
erage, not that the symmetry axes are somehow correlated 3 . 
Since there is no correlation in the signals between the two 
distinct lines-of-sight, there will be no contribution to the 
resulting CMB power spectrum from a clustering term. 



2.2 Two Sightlines 

For a single sightline through each outflow we found that 
the two-point correlation function vanishes, as the symme- 
try axes of LBGs are randomly oriented. When both lines- 
of-sight intersect the same LBG this cancellation does not 
take place. The average value of the two-point temperature 
anisotropy when both lines of sight intersect the same LBG 



AT V 



d 2 zP(z) 



d 2 h P(bi] 



d% P(b 2 ) ^ 
' (6) 



3 On small scales, the asphericity of the outflows may be corre- 
lated because they propagate into the same inhomgencous IGM, 
or because tidal gravitational forces produced correlations in the 
shapes of nearby galaxies iMackev et al.ll2002|y 



Performing the relevant integrations we find a non-zero an- 
swer when the product n e (r) v{f) has odd parity. The in- 
duced fluctuations produced when two lines of sight intersect 
the same LBG are due to Poisson fluctuations. As mentioned 
above, we parameterize the deviation from sphericity with a 
fudge factor e. Then the expectation value for two lines-of- 
sight intersecting the same LBG outflow region is 
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3 POISSON FLUCTUATIONS 

We have seen that there is an exact cancellation when the 
two lines-of-sight intersect two different LBG outflows and 
that a non-zero signal arises when the two lines-of-sight in- 
tersect a single LBG outflow. Correlations in the one LBG 
terms are produced by Poisson fluctuations in the number of 
intercepted LBGs. We will analyze this effect in two stages: 
first, we will consider the effect along a given line-of-sight 
(pencil-beam survey) and then generalize to the case of a 
finite beam size. 



3.1 Pencil Beam Survey 

The temperature anisotropics induced by LBGs of halo mass 
between M and M + dM, formed between redshifts Zf and 
Zf + dzf and scattering the CMB between redshifts z to 
z + dz is 
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where -Pjv(n LB G) i s the Poisson probability that n LB G LBGs 
are observed and dN is the mean number of LBGs in the 
redshift interval between z and z + dz, 
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(1 + z)H(z) dMdzf 

where R is the radius of an outflow at redshift z produced by 
an LBG of mass M formed at redshift Zf. The total signal is 
found by integrating over dM, dzf and dz. The expectation 
value toward a given line of sight vanishes since the expec- 
tation value for the temperature anisotropy produced by a 
single LBG vanishes. Nevertheless the variance of the sig- 
nal does not vanish due to Poisson fluctuations. For a single 
sightline through an LBG outflow region the contribution to 
the temperature anisotropies can be either positive or neg- 
ative with equal probability, however for two sightlines the 
contribution is always non-negative. The resultant variance, 
which is the case of two sightlines at a separation less than 
the characteristic angular size of LBG outflows, is 
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The rms value of the anisotropy can be evaluated in terms of 
the unknowns and the time changing Legendre coefficients 
ni and ve . For simplicity, we adopt the spherically symmetric 
solution for the shock radius and velocity, and parameterize 
the degree of asymmetry in the outflow by the fudge factor 
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3.2 Window Function Effects 

The average angular size of an LBG outflow #lbg ~ 10" 
is below the resolution of the upcoming generation of ex- 
periments, and so we must properly account for the beam's 
window function. The observed temperature anisotropy will 
be an average over a window function W(n) 
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For simplicity, we will take the window function shape to be 
a top hat of angular size 9, namely W(h) = 1/9 2 if jn| 9, 
and W(n) = if |n| > 9. 

The variance in an angular aperature defined by the 
window function, 
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In order for Poisson fluctuations to give a nonzero re- 
sult, the two lines of sight fi x and h' 2 must intersect the same 
LBG. Therefore they must be separated by less than #lbg- 
This requirement allows us to evaluate Eq. I12H as 
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Note that the Fourier multipole corresponding to #lbg is 
£ = 2-7t/#lbg- Here Da(z) is the angular diameter distance 
to redshift z. 

A broad window function allows for a larger number 
of LBGs within the beam. The window function is normal- 
ized such that it integrates to unity, and so one is observing 
the fractional fluctuations in the signal as an average over 
# 2 /$lbg independent coherence patches in the beam. This 
is equivalent to a fractional fluctuation of 1/ y/N as expected 
from Poisson fluctuations. 



4 RESULTS 

We numerically solve for the evoluti on of the shock front 
in the thin shell app roximation llTegmark et alJ Il993t 
iFurlanetto fe Loebl l2003'l . taking into account the effects of 
the halo gravity and the self gravity of the mass shell, the 
internal pressure of shocked IGM and the acceleration due 
to the cosmological constant. When calculating the inter- 
nal pressure we include Compton cooling, adiabatic cooling, 
and the addition of shock heated gas. Initially we assume 
that a fraction /„ =0.1 of the baryons assembled into the 
central LBG form stars with a Scalo initial mass function. 
In this c ase there is one supernov a per 126 Mq of star for- 
mation iFurlanetto fc Loebi 12003). Supernovae characteris- 
tically produce 10 51 ergs of energy, but only a small fraction 
of this energy, /sn, couples to the outflow, with the rest 
being radiated away. We adopt a value of /sn = 0.01; see 
IFurlanetto fc Loebl (2003) and reference within for a more 
detailed description of our outflow model. 



Assuming a Sheth-Tormen mass function 
JSheth fc Tormenl I1999T) we include the effects from all 
possible halos above the minimum galaxy mass. The lowest 
galaxy mass is determined by the maximum between the 
Jeans filtering mass and the cooling mass (dictated by 
the halo's ability to cool through atomic hydrogen line 
emission) 4 . We allow the LBGs to form between 5 z = 2 and 
z = 8. We continue to allow the outflow to evolve until its 
velocity equals to Hubble flow at its radius from the LBG or 
until z = 0. Increasing the upper redshift has little effect on 
our results because of the higher minimum LBG mass, as 
well as, the higher velocity of the Hubble flow which causes 
the LBG outflows to merge with the IGM at an relatively 
earlier time. Decreasing the minimum LBG mass could 
have a significant effect on our results because the ratio of 
the initial outflow velocity to the halo escape velocity at 
,the initial radius scales as v- m it/v e3c oc M~ 2 ^ . However, 
in low mass haloes the star formation timescale increases 
and supernova feedback becomes capable of decreasing /*. 
Some starburst activity continues to lower redshifts, albeit 
with a reduced intensity compared to high redshift star 
formation, changing the lower redshift of from zf = 2 to 
Zf = 1 would change our results by a factor of 2. Since the 
observed star formation efficiency decreases with decreasing 
zf, our model, which assumes that star formation and the 
subsequent supernova feedback only depends on the LBG 
mass, overestimates this change. 

In Fig. |5]we show the CMB power spectrum produced 
by several secondary mechanisms. The three blue dot-dashed 
curves denote the LBG outflow signal calculated in this Let- 
ter for the values of e = 1, 0.5, 0.25 from top to bottom. The 
green dotted curve delineates the Ostriker-Vishniac effect 
and the dashed cyan cur ve describes the patch y reioniza- 
tion effect as calculated bv lMcQ umn ct al. (2006) 6 Also, for 
comparison, we show the primary CMB anisotropics (solid, 
black curve) and the small scale C osmic Background Ima ger 
(CBI) data points (red triangles) llReadhead et alJliooH) . 

To characterize our results, let us mention some typical 
quantities for a common LBG halo of mass 5x 1O 9 M0 formed 
at zf = 2. In this case the outflow reaches a maximum 
comoving radius of 150 kpc at z = 0.9 before it merges 
with the Hubble flow. At late times the outflow velocity with 
respect to the local Hubble flow is (v — HR)/c ~ 7 x 10 -4 
and the Thomson scattering optical depth through the shock 
front isrsi 10" 5 . 

These small characteristic values for the Thomson scat- 
tering optical depth and outflow velocity imply that the sig- 
nal should not be notably polarized. Thomson scattering 
of a radiation field containing a quadrupole moment will 
produce polarized radiation. There are two standard ways 



4 Note that we include haloes with masses below the observa- 
tional sensitivity for LBGs, as there is no fundamental reason to 
exclude these haloes. 

5 We include this lower redshift limit because the phenomenon of 
downsizing, as well as the evolution of AGN luminosity function, 
implies that massive galaxies finished forming stars around that 
redshift. ^^^^^^^^^^ 

6 Note that iMcQ uinn et ah! l2006l> used a different cosmological 
model with a higher value of ag = 0.9. This artificially raises the 
amplitude of the results compared to value of ag = 0.73 adopted 
in this work in line with WMAP3. 
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Figure 2. Power spectra produced by several secondary mecha- 
nisms. The contribution from the LBG outflow is shown as the 
three blue dot-dashed curves for e = 1, 0.5, 0.25 from top to bot- 
tom. The primary CMB anisotropics are shown in the solid-black 
curve; the Ostriker- Vishniac and patchy r eionization contribu- 
tions as calculated by McQuinn et al. ( 2006) are dotted-green and 
dashed-light-blue respectively. The observed CBI data points are 
plotted as red triangles. 

in which scattering by a halo can produce polarization - 
(i) photons will double scatter in the halo; (ii) the peculiar 
velocity of the scatter er will induce a quadrupole moment 
in the radiation field dZel'dovich fc Sunvaevlll980tl . In the 
first case, the radiation can scatter in the halo producing an 
anisotropic radiation field; a second scattering of that radia- 
tion field can produce polarization at the level 0(t 2 v/c). In 
the second case, the peculiar velocity of the scatterer per- 
pendicular to its line-of-sight with respect to the observer 
will induce a quadrupole in the incident radiation field at 
the order 0(v 2 /c 2 ). A fraction r of the radiation field will 
scatter and become polarized at the level 0(rv 2 /(?)■ Since 
values of r and v/c are so small, we conclude that the po- 
larization power spectrum (which is sixth order in the small 
parameters of r and v/c) is negligible. 



5 DISCUSSION 

The contribution of outflows from Lyman-break galaxies to 
the CMB power-spectrum on arcminute scales is propor- 
tional to the square of their characteristic level of deviation 
from sphericity, e 2 . Future CMB experiments could therefore 
calibrate the intricate feedback process of galactic outflows 
on the IGM. Most tools used to study these feedback pro- 
cesses focus on the inner few kpc of the LBGs, even though 
the shock front is typically located at several of tens or hun- 
dreds of kpc. The secondary CMB anistropies calculated in 
this work provide a unique probe of these extended per- 
turbed regions around starburst galaxies at high redshifts. 

Even though the amplitude of the power spectrum pro- 
duced by this effect is small, the signal has distinctive spec- 



tral and spatial characteristics. The power spectrum has 
the same scaling with Fourier multipole (oc I 2 ) as radio 
or IR point sources. However, the frequency dependence of 
the anisotropies produced by our effect has the standard 
blackbody spectrum, whereas the radio point sources have 
a power law frequency spectrum. This distinct feature of 
our effect may allow future small scale experiments, many 
of which have excellent frequency coverage, to separate the 
anisotropies produced by LBGs from radio point sources. 
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